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Starting from HFB-6, we have constructed a new mass table, referred to as HFB-8, including 
all the 9200 nuclei lying between the two drip lines over the range of Z and TV > 8 and Z < 120. 
It differs from HFB-6 in that the wave function is projected on the exact particle number. Like 
HFB-6, the isoscalar effective mass M* is constrained to the value 0.80M and the pairing is density 
independent. The rms errors of the mass-data fit is 0.635 MeV, i.e. better than almost all our 
previous HFB mass formulas. The extrapolations of this new mass formula out to the drip lines do 
not differ significantly from the previous HFB-6 mass formula. 
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I. INTRODUCTION 

In the last few years we have been able to construct 
complete mass tables by the Hartree-Fock-Bogoliubov 
(HFB) method Q, SH H, with the parameters of the 
underlying forces being fitted to essentially all of the 
available mass data. This paper is the third in a se- 
ries of studies of possible refinements and modifications 
to our HFB- 2 mass formula PS, the first of our models 
that was able to give a satisfactory fit to the new data 
that had accumulated since the 1995 Atomic Mass Eval- 
uation 0. The most obvious reason for making such 
changes would be to improve the data fit, but there is 
also a considerable astrophysical interest in being able 
to generate different mass formulas even if no significant 
improvement in the data fit is obtained. The main point 
here is that the r-process of stellar nucleosynthesis pro- 
ceeds through the formation of nuclei that are so highly 
neutron rich that their properties cannot be measured 
but must be inferred by extrapolating the properties of 
known nuclei, and there is no guarantee that mass formu- 
las giving equivalent mass-data fits will still give the same 
masses out towards the neutron drip line. Moreover, even 
if they do, it is still possible that the underlying model 
(forces) will give different results for other properties rel- 
evant to the r-process (see Section 1 of the first paper of 
this series 0, referred to here as paper I). Predictions 
for the equation of state of neutron star matter could 
likewise differ. 

In paper I we discussed the role of a possible den- 
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sity dependence of the pairing interaction, while in pa- 
per II |j] we examined the question of the effective nu- 
cleon mass to be imposed on the Skyrme force. In the 
present paper we extend our HFB model by the restora- 
tion of the particle-number symmetry. In Section ^ we 
will recall the main features of the earlier HFB models 
relevant for the present study. In Sect. IIIII we present 
our formalism for particle-number restoration. In Sec- 
tion ^3 we describe a new fit of the force parameters to 
the mass data, with the resulting parameter set being 
labelled BSk8, and discuss its ability to predict masses 
and radii. A complete new mass table, HFB-8, is con- 
structed and compared with the earlier HFB-6 mass table 
in Sect. El Conclusions are drawn in Sect. IVll 



II. SUMMARY OF THE MODEL 

For convenience we recall some of the essential features 
of these HFB models. They are are based on conventional 
Skyrme forces of the form 

< = toil + xoPa)S{r,,) 

+ h{l + xiP^)^{pl5{r,,) + h.c.} 

+ t2(l + X2Pa)^Pi]-S{rij)pij 

+ ^Wo{(Ti + (Tj).p,j X S{rij)p,j (1) 

in the particle-hole (ph) channel, and in the particle- 
particle (pp) channel a (5-function pairing force acting be- 
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tween like nucleons treated in the full Bogoliubov frame- 
work 
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where p = p(r) is the local density, and po is its equilib- 
rium value in symmetric infinite nuclear matter (INM). 
Actually, it was only with models HFB-3 HFB-5 4j, 
and HFB-7 Q that the possibility of a density depen- 
dence in the pairing force was admitted; in all our other 
HFB models, HFB-1 0, HFB-2 1^, HFB-4 PJ, and HFB- 
6 we had ?7 = 0, as will be the case in the present 
paper. 

An important aspect relating to 5-function pairing 
forces concerns the cutoff to be applied to the space of 
single-particle (s.p.) states over which the force is allowed 
to act: both BCS and Bogoliubov calculations diverge if 
this space is not truncated |EI3- However, making such a 
cutoff is not simply a computational device but is rather a 
vital part of the physics, pairing being essentially a finite- 
range phenomenon. To represent such an interaction by 
a ^-function force is thus legitimate only to the extent 
that all high-lying excitations are suppressed, although 
how exactly the truncation of the pairing space should 
be made will depend on the precise nature of the real, 
finite-range pairing force. It was precisely our ignorance 
on this latter point that allowed us in Ref. to exploit 
the cutoff as a new degree of freedom: we found there 
an optimal mass fit with the spectrum of s.p. states £i 
confined to lie in the range 



Ef- £K< £t < Ef+£a, 



(3) 



where Ep is the Fermi energy of the nucleus in question, 
and £a is a free parameter. We shall adopt the same 
parametrization in the present paper. 

The mean-field states for even and odd particle num- 
ber have different structure. While they are HFB states 
for even particle number, they are one-quasiparticle ex- 
citations for odd particle number 



|MF) = |HFB) 
IMF) = aJlHFB) 



even particle number, (4) 
odd particle number. (5) 



Throughout this paper, we use the canonical basis to rep- 
resent the HFB states. The blocking of the unpaired par- 
ticle is not calculated completely self-consistently. Doing 
so requires breaking the time-reversal and axial symme- 
tries, and adding extra fields to the Hamiltonian. This 
is still too costly for a large-scale mass fit as performed 
here. Therefore we use an approximation, where the pair- 
ing correlations in odd nuclei are determined by removing 
the last nucleon 0. The occupation probability of the 
last occupied level is later corrected by adding the oc- 
cupation contribution of the unpaired nucleon. The last 
nucleon therefore does not contribute to the pairing (it 
does through the pairing tensor in the spherical config- 
uration), and the corresponding level should be treated 



accordingly. For the unprojected mean-field states, the 
occupation probability of the last occupied level k is 



(d(fc) - 2)v 



d{k) 



(6) 



where the degeneracy d{k) equals 2 for deformed nuclei 
and 2j -|- 1 for spherical ones (j being the angular momen- 
tum quantum number); vf, is the occupation probability 
obtained from the HFB equations for the level k without 
the blocked nucleon. To compensate for the approxima- 
tions done for odd particle number, we incorporate what 
is missing into the effective interaction and shall allow the 
pairing-strength parameter V-j^q to be different for an odd 
particle number (V^^) than for an even particle number 
(y^). For example, the pairing force between neutrons 
depends on whether TV is even or odd. We also allow 
the pairing-strength V-j^q to be different for protons and 
neutrons, as in all our previous mass models. Another 
feature of the HFB-2 model that we retain here is to 
add to the energy calculated by the HFB model a phe- 
nomenological Wigner term of the form 



Ew = Vw exp ■ 




(7) 



The method used to solve the HFB equations for Skyrme 
interactions has been presented earlier in Ref. 



III. RESTORATION OF THE 
PARTICLE-NUMBER SYMMETRY 

Mean-field approaches, such as the HFB used here, es- 
tablish an intrinsic frame of the nucleus and consequently 
break several symmetries of the Hamiltonian and the 
wave function in the laboratory frame ^ JJ. In par- 
ticular, finite nuclei break translational invariance, de- 
formed nuclei rotational invariance, reflection asymmet- 
ric shapes the parity symmetry, and the HFB framework 
the particle-number symmetry. These symmetry break- 
ings are required to include the desired correlations to the 
modeling (as multi-particle- multi-hole states), but at the 
same time gives rise to an admixture of excited states to 
the calculated ground state. The broken symmetries can 
be restored rigorously by projecting the wave function on 
the exact quantum numbers. A simpler procedure aims 
at estimating the contribution to the binding energy in a 
suitable approximation, and to add the resulting correc- 
tion to the binding energy. We adopted such a procedure 
in some of our previous mass formulas, in particular to es- 
timate the center-of-mass (cm) correction from the recoil 
energy, and the rotational correction within the crank- 
ing model [13. As far as the cm correction is concerned, 
the approximate prescription of Butler et al. [Tll | was 
replaced in Q by a more fundamental calculation of the 
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the recoil energy. However, so far no attempt to restore 
the particle-number symmetry has been undertaken in 
any of our HFB calculations related to global mass fit. 
The present paper is devoted to the projection of the 
wave function on the exact number of particles after a 
variation that includes the approximate Lipkin-Nogami 
projection before variation (referred to as PLN), and its 
impact on the mass fit. 

The pair correlations included in the mean-field HFB 
wave function |MF) are known to break the particle- 
number symmetry, as an HFB state is not an eigenstate 
of the particle- number operator. While the expectation 
value of the particle number (MF|iV|MF) is constrained 
to have the desired number N of nucleons on the aver- 
age via a Lagrange multiplier (the chemical potential), 
its dispersion 

(MF|(A7V)2|MF) = (MFjiV^lMF) - (MF|/V|MF)2 
= 2j2ulvl (8) 



Mr, 



11 for both protons and neutrons in all calcula- 



never vanishes in the presence of pairing. We use the 

\ is the occupation 
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standard notation, where vf, 
probability of the level k. 

The particle-number symmetry is restored by applying 
a particle-number projection operator on the mean- 
field state H in El El 



1$) 



P^|MF) 



(MF|P^|MF)i/2 



(9) 



which eliminates all components of the HFB state with 
a particle number different from N . |<1>) is a normalized 
eigenstate of the particle-number operator. The particle- 
number projector is given by 



pN ^ J_ 
2tt 



2tt 



^i4,{N-N) 



(10) 



The integration interval of the gauge angle (j} can be re- 
duced to [0, tt] for intrinsic wave functions with a well- 
defined "number parity" quantum number, i.e., for even- 
even systems as well as odd systems if treated in the 
blocking approximation (see below). For practical appli- 
cations, the integral over the gauge angle in the particle- 
number projection operator has to be discretized. 



P„(M,) = ^^e'^'"(*-"). 



(11) 



We use the Fomenko prescription 'l 5| , where 0m = 
7r(m — 1)/Mg and Mg is the number of angles in the 
interval [0, tt] used in the calculation. It can be shown 
that this particular choice of the projector removes all 
unwanted contaminations of the states up to n±2Mg j^. 
In order to avoid certain numerical problems that arise 
at (/) = f when has accidentally the value 0.5 0| we 
restrict ourselves to odd values of Mg, choosing finally 



tions. The particle number dispersion is subsequently 
checked to verify that the projection actually worked. 

The particle-number symmetry restoration is per- 
formed in two steps. First, an approximate projection 
before variation is performed within the Lipkin-Nogami 
prescription ES E^ E^ EM E3 1 where we take only the 
pairing part of the interaction into account when calcu- 
lating A2. In a second step, the converged HFB states 
are then exactly projected as described by Eq. ifTUI) . As 
usually done, we constrain the HFB states to the same 
particle number that we project out exactly afterwards. 
In such a projection after variation approach, the main 
task for the Lipkin-Nogami scheme is to enforce the pres- 
ence of pair correlations in the HFB states for all 
and Z and all deformations to avoid the breakdown of 
pairing for small level densities. This two-step approach 
was shown to be an excellent approximation to the ex- 
act self-consistent particle-number projection before vari- 
ation Ell- 

The nuclear mass and all other observables are cal- 
culated from the projected state applying a generalized 
Wick's theorem. The basic contractions needed for that 
are given by E3 



Pkk'i<P) = 



Kkk' {<P) 



i^kk' (0) 



(MF|4afc'e*"*^|MF) 



(MF|MF((/.)) 
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^kk' = Pkk'W , 



(MF|44,e*'^*|MF) 
(MF I MF (</>)) 

-UkVk ^ 

(MF|afcafc,e^'^^|MF) 
(MF|MF(</>)) 



UkVke 
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Skk' = -Kfcfc' ((/>), 



(12) 



(13) 
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where we introduced the shorthand notation |MF((/))) = 
e"^^|MF) for a MF state rotated in gauge space. We 
employ the phase convention {uj.,vj.) — {uk, — Wfc), where 
k denotes single-particle states with positive angular mo- 
mentum projection on the symmetry axis, and k the 
corresponding time-reversed states with negative angu- 
lar momentum projection on the symmetry axis. For the 
blocked orbital in a mean-field state with odd particle 
number, the contractions of the level k are determined by 
removing the blocked nucleon, since it should not enter 
the particle number projection. Of course, the contribu- 
tion of this unpaired nucleon is added to Eq. H12I) after- 
wards (see, for instance, Eq. (|23l) ). In order to express the 
expectation value of observables in the projected state, 
it is useful to introduce normalized 0-dependent overlaps 



O(0,n) 



(MF|MF(0)) 
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E™ti(MF|MF(0„O)e-™'^" 



(15) 
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Based on Eq. (|12(l . the number of particles of the pro- 
jected state can then be written as 



($|7V|$) = ^O(0™,n)Tr[p(0,„ 



which reduces to the famihar {N) = 2 ^ 



(16) 



The particle-number projected energy depends on the 
0-dependent generalized densities and is calculated by an 
integration over the spatial coordinates r and the sum- 
mation over the two sets of gauge angles (j)'^ {q — n,p). 
The total mean-field energy can therefore be expressed 
as 



E„ 



rriji — l m-n — l 



(17) 



where Smf is the Skyrme energy density functional (ob- 
tained with Eq. (3a) of Jl^] and the (/>— dependent densi- 
ties). The pairing energy is calculated similarly: 



Ma 



E 

rnn— 1 mp — 1 



E 



where 



pair,g 



(r, 



X / dr £pair,g(r,C„:CJ: 



= ^,,Tr[K(,/)«)]Tr[K 

"(/j(r, f/-: 



(18) 
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For the density-dependence in E^^^f and -Epair 
the mixed density obtained from Eq. (|12|l . as done in 
most applications involving the mixing of different mean- 
field states. This choice, however, is not unique, see the 
discusion in Refs. |13Ll20| and references given therein. 



1. Decomposition of mean-field states 

The weight of the n'-nucleon wave function in the 
decomposition of the unprojected n-particle HFB state 
|MF„) in given by the overlap 



On' 



(MF„|P" IMF, 



(20) 



This quantity can be used to illustrate the particle num- 
ber dispersion of the HFB states. 

An example is given in Fig. ^ for ^''"Pu, where we 
compare the spherical configuration with the deformed 
equilibrium one, which has a quadrupole moment of 
Q = 2782 fm^. The decomposition is significantly dif- 
ferent. In lowest order, the distribution of the weights 
has a Gaussian shape 



((MF|iV|MF) -n'f 
2^ 



On' 



V27rcr2 



exp 



, (21) 




140 150 160 
Z N 
FIG. 1: Proton (left) and neutron (right) overlaps between 
the unprojected HFB wave function of "^""^Pu and its state 
projected on n {N or Z) nucleons; diamonds refer to spherical 
overlaps, crosses to deformed ones; the full line corresponds 
to a Gaussian centered on n = 94 (protons) or 146 (neutrons) 
and of standard deviation equal to the particle-number 
dispersion. 



where the width of the Gaussian cr^ = (MF|(AiV)2|MF) 
is determined by the particle-number dispersion of the 
unprojected state, Eq. ©. In the specific case of ^''"Pu, 
the Gaussian is more spread for the spherical shape, 
where pairing is stronger due to the larger level density 
around the Fermi surface. 

The different blocking prescription in spherical and 
deformed configurations leads to significant variations 
in the weight of neighbouring nuclei in the wave func- 
tion, as seen in Fig. |21 for the spherical ^^In and ^"^"in 
isotopes: the blocked nucleon in the spherical case is 
not removed from a doubly-degenerate level (to leave it 
empty) but removed from the gg/2 level, which is equiv- 
alent to the removal of 20% of a nucleon on the last 
five doubly-degenerate levels of a deformed configuration. 
The weight of the wave functions of adjacent nuclei is 
therefore higher (in the Z = 49 or iV = 51 In isotopes) 
if the spherical code is used. It is also observed that the 
deviation between Eq. H21|l and Eq. (|20|l is more pro- 
nounced, due to a shell closure effect. 

For consistency, the cm correction should be calcu- 
lated from the particle-number projected state. Written 
concisely, the gauge angle-dependent cm energy of the 
particle-number projected state reads 

£cm{(i)) = [2 E Pfefc Afcfc + 2 E (PkkPll + Kkkiill) 



k>0 



k,l>0 



For odd nuclei, {pkkPii + i^kki^u) becomes 
-{d{k)-25b{k))pkk + 5b{k) 



(22) 



d{k) 

{d{l)-2St{l))pn+Sb{l) 
d{l) 



+ '^kk^U 



(23) 



with (5ft (fc) = 1 if the blocked nucleon is on the level k (0 
otherwize). The cm energy (for one nucleon species q) is 
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45 50 55 
N 

FIG. 2: Same as Fig.^for the spherical ^^In and ^""in istopes. 
The Gaussian approximation, Eq. 1211 (thin line with dots or 
dashes obtained from the spherical or deformed code) can- 
not describe the pronounced shell closure. For Z = 49 or 
N = 51, the overlaps obtained with the deformed code differ 
significantly from the overlaps calculated assuming spherical 
symmetry, due to a different implementation of the blocking 
approximation. 



then 



Mo 



(24) 



Finally, note that as far as the rotational correction is 
concerned, we will restrict ourselves to the simplified 
cranking formula, although ideally the rotational sym- 
metry should also be restored with projection techniques. 
So far, in all our global mass fits, we adopted the rota- 
tional correction calculated from the HFB states 



(MF|J2|MF) 
21 



(25) 



where J is the total angular momentum operator in the 
intrinsic frame and T the moment of inertia. The crank- 
ing model gives this latter quantity as 



X.. = 2 E 



\J.\k') 



k,k'>0 



Ek + Ek' 



-{ukVk> - Uk'Vkf, (26) 



where the summation runs over quasi-particle s.p. states, 
the Ek are the corresponding quasi-particle energies, the 
matrix elements are calculated in the canonical basis, and 
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are the corresponding occupation probabili- 



While the cranking model is in general quite satisfac- 
tory, a major problem in our mass-model applications 
arises from the fact that a priori we do not know which 
nuclei will be spherical (or quasi-spherical), for which nu- 
clei a numerical problem arises from the 0/0 indetermi- 
nacy in Eq. H25|l : Erot must vanish, of course, for spher- 
ical nuclei. A way around this problem, introduced in 
Ref. [2JI , is to take for the X the linear combination 



aXcr + (1 — a)X, 



rtg 



where Xrig is the rigid-body moment of inertia 



M 



1 



-i2R'A+-Q) 



(27) 



(28) 



in which R is the rms matter radius of the nucleus, M 
the nucleon mass, and Q the quadrupole moment. The 
value of a was originally taken to be 0.8, but was reduced 
to 0.75 in Ref. [23 and kept unchanged in all subsequent 
mass formulas (HFB-l-HFB-7). However, we see from 
Tabled that with this prescription the deformation ener- 
gies AEiso of highly deformed shape isomers, i.e., their 
energy relative to the ground state of the nucleus in ques- 
tion, are badly underestimated, and can even be negative. 
For this reason we adopt for the rotational energy the 
phenomenological prescription 



Er 



^j^crank 



" tanh(c|/32|) , 



(29) 



or equivalently introduce a new prescription based on 
the cranking value of the moment of inertia modified as 
follows: 



rth = ■^Xc.rCOth(c|/32|). 



(30) 



ties. 



Here, the dimensionless quadrupole moment (32 is defined 
as a function of the quadrupole moment Q2 and the re- 
duced radius Rq by (32 = \/5TrQ2/{3ARQ). Experimental 
information on the mass of well deformed nuclei, but also 
the energy of the shape isomers in the actinide region is 
used to determine the free parameters (&, c). Values of 
(&,c) range from (0.6,5.5) to (0.7,3.5). For HFB-8, we 
adopt b = 0.65 and c = 4.5. We find that this pre- 
scription not only is equivalent to the "mixed" prescrip- 
tion H27|l for ground states, i.e., for masses, and likewise 
avoids any problems in the spherical limit, but also leads 
to considerable improvements in the estimates for the 
shape isomers, as can be seen from the penultimate col- 
umn of Table n 

The rotational correction should be accompanied by 
corrections for vibrational zero-point motion. To the best 
of our knowledge, there exists no strategy to estimate 
them properly and at the same time include them in a 
global mass fit as ours with current computing resources. 
Common strategies are to calculate RPA correlations [s^, 
Isij . usually restricted to spherical nuclei, or to use the 
generator coordinate method jSJ. More research in this 
direction is necessary in the future. 
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TABLE I; Comparison between experimental I2 I I27 I I2I . I2I 

and theoretical (reflection symmetric assumed) energies of 
shape isomers obtained with the HFB+PLN(BSk8) model 
using three different rotational correction prescriptions: Eq. 
lOni (modified cranking), Eq. (1771 (mix with a = 0.87), and 
Eq. 1281 (rigid). The last two lines show the rms and mean 
deviations, respectively, compared to experiment. All quan- 
tities in MeV. 



Z N A cxp crth mix rig 

2.2 
2.1 
2.3 
1.9 
1.8 
1.9 
2.0 
1.5 
1.3 
1.4 
1.1 
1.4 
1.3 
1.3 
1.5 
0.9 
1.6 
1.7 
1.0 
1.1 
1.3 
0.6 
1.3 
1.3 
0.7 
0.8 
0.6 
1.0 
0.9 
0.1 



o- 0.64 1.62 0.95 

e 0.40 1.54 0.80 



IV. THE MASS FIT 

To study the impact of the PLN framework, we start 
from the BSk6 force P|, keeping some of its character- 
istics: an isoscalar effective mass M*/M constrained to 
0.8, the symmetry energy J to 28 MeV, the 7 exponent in 
the ^3 term of Eq.(^) to 1/4, and a density-independent 
pairing interaction (i.e., rj — 0). The results of the HFB- 
8 mass fit of the new HFB-l-PLN model are presented in 
Table ^1 where for comparison we also show the results 
of the HFB-6 model. Since the latter was fitted to the 
2135 nuclei with Z,N > 8 whose masses had been mea- 
sured and compiled in the unpublished 2001 Atomic Mass 
Evaluation (AME) of Audi and Wapstra 33] wc show re- 
sults for both this data set and the 2149 measured masses 
of the updated AME that was released and published at 
the end of 2003 [23| ■ Experimental and calculated mass 
excesses are compared in Fig. 

The parameters of the force BSk6 and BSk8 are com- 
pared in Table imi while the macroscopic parameters, i.e., 
the parameters relating to (semi-) infinite nuclear matter 
calculated for these forces, are shown in Table HVl These 



TABLE II: Rms (a) and mean (e) errors (in MeV) in the 
predictions of masses M obtained with the BSk6 and BSk8 
forces with respect to the 2135 nuclei in the 2001 AME [s^ 
and the 2149 nuclei in the 2003 AME 0. The last two 
lines correspond to the rms and mean errors (in fm) in the 
predictions of the 523 measured charge radii (tc)- 





BSk6 


BSk8 


a{M) (2135 nuclei) 


0.684 


0.659 


e(M) (2135 nuclei) 


0.021 


0.005 


a{M) (2149 nuclei) 


0.666 


0.635 


e(M) (2149 nuclei) 


0.014 


0.009 


o-{rc) (523 nuc) 


0.0262 


0.0250 


e(rc) (523 nuc) 


-0.0028 


0.0047 



4 

3 



2 - " 




M(Exp)-M(HFB8) 

.4 f , , ^ I - 

20 40 60 80 100 120 140 160 

N 

FIG. 3: differences between experimental and calculated 
mass excesses as a function of the neutron number for the 
IIFB-8 mass tables 



include, in addition to the previously defined parameters, 
ay, the energy per nucleon at equilibrium in symmetric 
infinite nuclear matter, M*/M the ratio of the isovector 
effective nucleon mass at density po to the real nucleon 
mass M, K^, the incompressibility. Go and Gq the Lan- 
dau parameters as defined in Ref . [35| , p frmg , the density 
at which neutron matter flips over into a ferromagnetic 
state that has no energy minimum and would collapse 
indeflnitely 36], Usf, the surface coefficient, and Q, the 
surface-stiffness coefficient |0] (the meaning of Go, Gq, 
and p frmg is critically discussed in ) • 

The Skyrme parameters have changed little from the 
original BSk6 Skyrme force; the pairing strengths have 
instead decreased significantly, by 2% for T4„ and 11% 
for l/jrp. For the first time we obtain V^rp < V-rrm which is 
what one expects when the zero-range pairing force with 
different coupling constants for protons and neutrons. 
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1.04 


245 


1.7 


244 


2.0 



2.6 


1.6 


2.4 


1.3 


2.5 


1.5 


2.3 


1.2 


2.3 


1.2 


2.2 


0.9 


2.4 


1.4 


1.9 


0.6 


1.7 


0.7 


1.6 


0.6 


1.5 


0.4 


1.9 


0.8 


1.8 


0.7 


1.8 


0.6 


1.9 


0.8 


1.3 


0.1 


2.0 


0.8 


2.2 


1.0 


1.4 


0.4 


1.4 


0.4 


1.7 


0.6 


1.0 


-0.2 


1.7 


0.6 


1.6 


0.5 


1.1 


0.0 


1.3 


0.2 


1.1 


-0.2 


1.5 


0.3 


1.3 


0.0 


0.6 


-1.0 
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TABLE III: Skyrme-force and pairing-force parameters of 
BSk6 and BSk8 





BSk6 


BSk8 


to [McV fm ] 


ori/i Q Q 1 T 
-zU4o.oi / 


-zUoo.ozo 


tl [iViev im J 




oyo.ozuo 


62 [MeV im J 


1 70 Q7Q 

- i ( 0.0 / y 


-lyD.uuoz 


63 [iVlc V rm J 


1 9Pi1 1 7 


IZ-rtOO-OO 


Xq 


U. 1 OuOUc/ 


77'^^^9^^ 

U. 1 i OOiiO 


Xi 


n 7QQ1 




•1^2 


\J . tj L> L> 




X3 






Wn fA/ToV fm^l 
K r (J 1 iViC V Llll J 


142.4 


147.8 


! 


1 /4 


1 /4 


[ivie V im J 


-oz 1 .z 


-014:. U 


[iVio V rm J 


-00 i . y 


^9Q ft 


Vp [MoV fm ] 


on A c 
-oz4.0 


-293.0 


Vp [JVlcV rm J 


Q/10 /I 
-04Z.4 


-ouy.y 


V 








a 








£A [MoV] 


17 


17 


Vw [McV] 


1.76 


1.85 


A 


700 


780 


[MoV] 


0.58 


0.66 


Ao 


28 


26 


a (Eq.lEl) 


0.75 




6,c (Eq.lHU) 




0.65,4.5 


TABLE IV: Macroscopic parameters of the forces BSk6 and 


BSk8 








BSk6 


BSk8 


a„ [MoV] 


-15.749 


-15.824 


PO [fm-3] 


0.1575 


0.1589 


J [MeV] 


28.0 


28.0 


m;/m 


0.80 


0.80 


M„7M 


0.86 


0.87 


ifi, [MoV] 


229 


230 


Go 


0.066 


0.042 


Go 

Pfrmg/pO 


0.312 


0.261 


1.82 


1.70 


Usf [MeV] 


17.18 


17.64 


Q [MoV] 


53.4 


53.9 



Eq. lO, shall simulate an attractive isospin- invariant 
isovector pairing force and a repulsive isospin symmetry- 
breaking pairing force originating from the Coulomb in- 
teraction. The new Skyrme parametrization globally im- 
proves the prediction of experimental masses, as shown in 
Table im The comparison between the HFB-6 and HFB- 
8 mass tables, in particular towards the neutron dripline 
is further discussed in Sect. Ivl 

Table lll| also shows the rms and mean deviations be- 
tween theoretical and experimental charge radii for the 
523 nuclei listed in the 1994 compilation ^| (for more 
details on the HFB derivation of the charge radii, see 




r [fm] EXP 

ch 



FIG. 4: Comparison of the theoretical and experimental 
charge radii for the 523 nuclei listed in the 1994 compilation 
of H 




FIG. 5: Comparison of the theoretical (full line) and exper- 
imental (crosses with error bars) charge densities for for ^^S 
and ^'^*Pb. Experimental data are from |40|. 



|4ll|). The comparison is given in Fig. ^ The overall 
agreement with experiment is seen to be excellent, al- 
though some fine structure related to local anomalies is 
not reproduced in all details. Similarly, the radial charge 
density distribution of spherical nuclei as obtained from 
the charge form factor is found to be well reproduced 
even down to the centre of the nucleus, see the examples 
of ■^■^S and ^"^Pb given in Fig.jS] Such a good agreement 
between theoretical and experimental densities and radii 
is strongly related to the adopted value of the saturation 
density po (Table III) or equivalently the Fermi momen- 
tum kp = 1.33 fm~ given by (37r2po/2)i/3^ However, 
in the case of densities, the excellent agreement with ex- 
periment both in the surface and inside the nucleus must 
be regarded as fortuitous, being completely beyond our 
control in mass fits of the kind described here. 



8 



4 - 



> 



<1 -2 









4 




♦ *\*+ 


















M(HFB-8)-M(HFB-6) 
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FIG. 6: Differences between the HFB-8 and the HFB-6 
masses as a function of the neutron number for all nuclei 
with 8 < Z < 120 lying between the two drip lines. 
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FIG. 7: Differences between the HFB-8 and the HFB-6 
masses as a function of the neutron separation energy S„ for 
all nuclei with 8 < Z < 120 lying between the two drip lines. 



V. EXTRAPOLATIONS 

With the BSk8 Skyrme force determined as described 
we constructed a complete mass table, labeled HFB-8, for 



the same nuclei as were included in the HFB-6 tables, i.e., 
all the nuclei lying between the two drip lines over the 
range of Z and iV > 8 and Z < 120. 

The HFB-8 masses are compared in Figs. 16171 with the 
HFB-6 masses. 

Differences seldom exceed some 2 MeV, even close to 
the neutron dripline; the largest deviations are found for 
open shell nuclei. Shell effects far away from stability are 
found to be very similar for the two mass tables. More- 
over, the HFB-8 shell gaps are very similar to the shell 
gaps obtained with the HFB-6 mass formula so that they 
are not shown here. These results once again confirm 
the relative stability of the HFB mass predictions with 
respect to different parametrizations or frameworks, as 
already emphasized in 



VI. CONCLUSIONS 

The pairing correlations included in the HFB wave 
function are known to break the particle-number sym- 
metry. The restoration of the exact particle number is 
done on the basis of the projection technique, i.e by pro- 
jecting the wave function on the exact number of particles 
after a variation that includes the approximate LN pro- 
jection before variation. Doing so, we have constructed 
a new Skyrme force, labelled BSk8, the parameters of 
which reproduce the 2149 measured masses with an rms 
error of 0.635 MeV. The final table, referred to as HFB-8, 
includes all the 9200 nuclei lying between the two drip 
lines over the range of Z and N > 8 and Z < 120. The 
extrapolations of this new mass formula out to the drip 
lines do not differ significantly from the previous HFB- 
6 mass formula obtained without the restoration of the 
particle-number symmetry. 
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